Persistent currents in a graphene ring with 
armchair edges 

Bor-Luen Huang 

Department of Physics, National Taiwan Normal University, Taipei 11677, Taiwan 
Ming-Che Chang 

Department of Physics, National Taiwan Normal University, Taipei 11677, Taiwan 
Chung-Yu Mou 

Department of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan 
Institute of Physics, Academia Sinica, Nankang 11529, Taiwan 

Physics Division, National Center for Theoretical Sciences, P.O.Box 2-131, Hsinchu, 
Taiwan 

Abstract. A graphene nano-ribbon with armchair edges is known to have no edge 
state. However, if the nano-ribbon is in the quantum spin Hall state, then there must 
be helical edge states. By folding a graphene ribbon to a ring and threading it by a 
magnetic flux, wc study the persistent charge and spin currents in the tight-binding 
limit. It is found that, for a broad ribbon, the edge spin current approaches a finite 
value independent of the radius of the ring. For a narrow ribbon, inter-edge coupling 
between the edge states could open the Dirac gap and reduce the overall persistent 
currents. Furthermore, by enhancing the Rashba coupling, we find that the persistent 
spin current gradually reduces to zero at a critical value, beyond which the graphene 
is no longer a quantum spin Hall insulator. 

PACS numbers: 81.05.ue, 61.72.J-, 71.15.-m 
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1. Introduction 

A crucial feature of topological insulator is the existence of metallic surface states pQ. 
These surface states (or edge states if the insulator is two dimensional) have odd pairs 
of channels, in which each pair consists of opposite spins moving in opposite directions 
(helical edge states). Therefore, they could carry charge as well as spin currents. We are 
interested in studying their transport using a lattice model. In theory, one could rely on 
a model with external leads and reservoirs, or with a closed ring in which the currents 
are driven by a threading magnetic flux. In this work, we adopt the latter (simpler) 
approach, and use the graphene lattice as a theoretical model to study the transport of 
these important helical edge states. 

Graphene is one of the earliest candidates of two-dimensional topological insulator 
(also called quantum spin Hall insulator) [2J. It is found that, graphene would become 
a quantum spin Hall (QSH) insulator if there is a spin-orbit interaction (SOI). However, 
since graphene's SOI is only of the order of 10~ 3 meV [3], it poses an experimental 
challenge to probe such a state [3J. Experimental issues aside, graphene lattice is ideal 
for theoretical investigation of the helical edge states. To avoid complications from pre- 
existing edge states, we choose a graphene ribbon with armchair edges, which has no 
edge state in the absence of SOI. That is, the edge states studied here are intrinsic to 
the QSH insulator. In contrast, a graphene ribbon with zigzag edges has edge states 
even without SOI [5J. These edge states would become helical edge states (similar to 
those in this work) with SOI. 

If the graphene ribbon is wrapped up to form a tube (or a ring, with armchair 
edges), with a magnetic flux passing through the hole, then the helical edge state 
provides a robust channel for persistent charge current (PCC) and persistent spin current 
(PSC). PCC in a mesoscopic metal ring was first predicted [6], and later verified in 
experiments [7] decades ago. In a recent work [8], the PCC was measured accurately by 
mounting nano-rings on a micro-cantilever in an external magnetic field. In addition 
to PCC, the PSC in a ring in a textured magnetic field [9], or in a ring with SOI [10] 
has been proposed. Further researchers show that it is possible to have a PSC in a 
ferromagnetic ring [TT] or an antiferromagnetic ring [12] . 

There have been several studies of the Aharonov-Bohm effect and PCC in a 
graphene ring with different geometrical shapes (without spin-orbit interaction). For 
example, a nano-torus [13], flat rings (similar to Corbino-disks) with various types of 
boundaries [HI [151 [16], or a tube with zigzag edges [17]. Observation of the Aharonov- 
Bohm effect is an indirect evidence of the PCC. Recently, the Aharonov-Bohm oscillation 
in a nano-rod of topological insulator has been observed [18] and studied [191 ED]- 
Theoretical investigations of the persistent currents in a two-dimensional topological- 
insulator ring, which are based on an effective continuous model, can be found in 
Refs. [HI [22] . In Ref . [23] , the zigzag edge of a graphene sheet is found to be magnetized 
due to the Hubbard interaction. Since the time reversal symmetry is spontaneously 
broken, edge charge current appears without the need of a driving magnetic flux. 
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In this paper, we study the quantum spin-Hall states and the persistent currents 
of an armchair ring with SOI. Different from earlier approaches J2TJ [22], this work is 
based on a lattice model, which allows insulating band structure and genuine edge 
states. By changing its parameters, the gaphene ring can be in or out of the QSH 
phase, and the behavior of the persistent currents across the phase boundary can be 
studied. It is found that: For a broad ribbon, the edge PCC is small (but non-zero), 
while the edge PSC approaches a finite value. For a narrow ribbon, the edge states 
from opposite sides couple with each other, which may open an energy gap and reduce 
the persistent currents. We have also considered the Zeeman coupling and the Rashba 
coupling. The former introduces addition jumps in the periodic saw-tooth curves of the 
persistent currents. By breaking spin conservation, the Rashba coupling reduces the 
edge spin current, which eventually vanishes at the boundary of the QSH phase. Many 
of the results obtained here should also apply qualitatively to the helical edge states and 
persistent currents of a three-dimensional topological-insulator ring. 

The paper is outlined as follows: In Sec. 2, we describe the theoretical model being 
used. In Sec. 3 we first study the helical edge states and related PCC and PSC in 
rings with various sizes, then comment on the effects of the Zeeman and the Rashba 
interactions. Sec. 4 is the conclusion. 

2. Theoretical formulation 

We start with a tight-binding model for a graphene ribbon. One imposes the open 
boundary condition with armchair edges along the x-direction, and periodic boundary 
condition along the ^/-direction (see Fig. [1]). A magnetic flux $ passes through the inside 
of the ring along the (— x)-direction. Because of the magnetic flux, an electron circling 
the ribbon once would acquire an Aharonov-Bohm phase of 27r$/$ , where $o is the 
flux quantum. The Hamiltonian reads [2], 

H= -t e iei 4 Cj + i\ so Yl e^V^Vc, (1) 

<i,j> 

+ ^fi B B^2cla x Ci + iX R Y e iei c\((T x d^Cj, 

i <i,3> 

where cf = (cL, c|i ) creates an electron at lattice site i. The first term accounts for 
the nearest-neighbor hoppings. The second term is a SOI with next-nearest-neighbor 
hoppings, in which = (2/v^3)(di x d 2 ) ■ z, d\ and d 2 are the two nearest-neighbor 
bonds that connect site-j to site-i, and cr are the Pauli matrices. The third and the 
fourth terms are the Zeeman and the Rashba couplings. In the following, all energies 
and lengths will be in units of t and a (lattice constant) respectively. 

The phases 8\ and 8 2 are related to the magnetic flux: 8i = 20/3 for the nearest- 
neighbor bonds along the y-direction (see Fig. [T|), while 8i = <p/3 for the zigzag bonds, 
where = (27r/A 2/ )$/$o, N y is the number of unit cells around the ring. The phase 
82 is zero for the next-nearest-neighbor bonds along the x-direction, while 82 = <fi for 
the other next-nearest-neighbor bonds. These phases are chosen in such a way that the 
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a) 





Figure 1. The honeycomb lattice is shown in (a). Red and blue arrows are the 
next-nearest-neighbor hoppings. The red and blue arrows in the plot are for spin-up 
electrons in A- and i?-sublattices, and + and — signs refer to the signs of i/y. For 
spin-down electrons, red and blue arrows appear in opposite directions. The effective 
two-lag ladder for the Hamiltonian in Eq. (3) is shown in (b). 



electron would pick up the correct Aharonov-Bohm phase when circling the ring once, 
no matter which path it takes. The geometric curvature of the tube is known to enhance 
the SOI |24j, but such a curvature effect is not considered. 

Because of the translational invariance along the ^-direction, one can perform the 
following Fourier transformation 



,ik y y 



and map the system to a two-lag ladder 

as, 
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In Eq. (3), hbB is multiplied by an unit matrix. The Rashba term gives 
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-Me^'^ and = -^A^e 



i<t>k/3 



At zero temperature, the PCC of the nano-ring is given by (in units of et/h), 

I = > — —. (7) 



3$ 



in which one sums over filled states (for total current) or edge states (for edge current). 
The PSC for spin- z will be calculated by the following semiclassical expression [10] (in 
units of (h/2)t/h = t/Air) 



inky 



(8) 



where \mk y ) are the Bloch states. 

When the spin is conserved (Sees. 3.1 and 3.2), the following notations are used: 
The current for spin-up and spin-down electrons are I a . Therefore, the total charge 
current is / = I + + 7_, and the total spin current is I s = I + — The spin-dependent 
edge currents are I Ra and Iia- For example, the charge and spin currents for the right 
edge are I R = I R+ + I R ^ and I S R = I R+ - I R _ respectively. 
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Figure 2. Energy spectrum for a graphene ribbon with n x = 100, n y = 100, and 
Aso = 0.01. The inset shows the magnified view near the Dirac point. The bulk gap 
is about 0.136, which is slightly larger than the bulk value of 6\/3Aso because of the 
finite-size effet. The crossing edge states inside the gap is the main characteristic of 
the quantum spin Hall phase. 

In general, I a is the sum of edge current and bulk current. The existence of 
persistent bulk current in an insulator ring (the usual type or the QSH type) is due to 
the discreteness (due to the finite circumference) and asymmetry (due to the magnetic 
flux) of the energy spectrum. As a result, it only exists in a ring with finite radius (with 
or without SOI), and diminishes when the ring is larger. On the contrary, there is no 
spin-resolved current (either bulk or edge) in an usual insulator ring. It only exists in 
the QSH phase and is not reduced in a larger ring. That is, the persistent spin current 
reported below is intrinsic to the QSH phase. 

3. Numerical Results 

There are several tunable parameters in the numerical calculations: width and length 
(circumference) of the ribbon, spin-orbit couplings Xso and Xr, and the ^-factor (for the 
Zeeman coupling). In addition, the magnetic flux inside the ring controls the magnitude 
of the persistent currents. To explore such a large parameter space, the following path is 
adopted: we first show the energy spectrum and the decay length of the edge states, then 
present the charge and spin transport for broad ribbons (Sec. 3.1). Subsequently, the 
transport in narrow ribbons are shown and the finite-size effect discussed (Sec. 3.2). In 
either case, to simplify the discussion, Zeeman and Rashba couplings are not included. 
Their effect will be investigated in Sec. 3.3 and Sec. 3.4 respectively. Among these 
parameters, the most important ones are Xso and Xr, which dictate whether the 
graphene is in the QSH phase [2]. Their overall influence on edge states and spin 
current can be seen in Figs. 3 and 10. 

Before showing the result, we review the energy spectrum of a flat graphene nano- 
ribbon with armchair edges, but without SOI [5]. First, the width of the ribbon has a 
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Figure 3. Decay length of the edge state at the Dirac point. The size of the ring is 
n x — 100 and n y = 100. The decay length is inversely proportional to the strength of 
the spin-orbit interaction. 

major effect on the energy gap. If n x = 3n + 2, then there is a Dirac point at k y = 
and the system is gapless in the tight-binding model. If n x ^ 3n + 2, then the system 
has a finite gap, which approaches zero as n x becomes infinite. Second, in the presence 
of SOI, the Dirac point of a ribbon with n x = 3n + 2 is opened, and mid-gap edge states 
are crossing at k y = 0. Similar mid-gap edge states exist for n x ^ 3n + 2. 

3.1. Broad ribbon 

In Kane and Mele's analysis [2], the graphene enters the QSH state once the SOI is 
turned on, without a finite threshold value. Therefore, the edge modes for the QSH 
state should immediately appear, no matter how small the SOI is. In Fig. 121 we show 
the energy dispersion for Xso = 0.01 and $ = 0. Overall, the energy spectrum shows 
intangible change from the one with Xso — 0. However, in the inset one can see that 
the Dirac point is opened up, with mid-gap edge modes. The edge modes do not merge 
into bulk bands up to near the boundaries of the Brillouin zone. 

The envelope of the probability distribution for the edge state decays exponentially, 
^(x)! 2 oc exp(— d/£), where d is the distance from the edge. As one increases the SOI, 
the edge states should be more localized. In Fig.[3l we show the decay length £ of the edge 
state at the Dirac point for different values of Xso- Upon fitting, the inverse of the decay 
length is about 5.89A,so, proportional to the SOI. For example, for Xso = 0.01(0.05), 
the decay length is about 17 sites (3 sites). With this result, one can estimate whether, 
for a given ribbon, the inter-edge coupling could be safely neglected . 

At first, we report the PCC without SOI. In this case, the current is from bulk states 
only since there is no edge mode. The PCCs are shown on the left of Fig. HI As we 
mentioned earlier, there is a Dirac point in the energy spectrum when n x = 3n + 2. The 
discontinuity of the slopes at the Dirac point produces the finite jump at zero flux for 
n x = 50. On the other hand, systems with n x ^ 3n + 2 are gapped (without the Dirac 
point), so the PCCs show no jump for n x = 51, 52. Furthermore, the overall magnitude 
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Figure 4. The plots on the left column are the persistent charge currents for A50 = 0, 
and those on the right column are for Aso = 0.05. The charge currents are in units 
of et/h (see Eq|7]).The figures on the first, the second, and the third rows are for 
n x = 50, 51, and 52 respectively. On the right column, one can see that the persistent 
currents decrease as n y becomes larger. The magnitudes for the 3 plots on the right 
are roughly the same, in sharp contrast to the ones on the left. 



of the PCC is also much smaller. 

The plots for the PCCs with SOI are shown on the right of Fig. HI Even though 
the SOI introduces little change to the PCC for n x = 50, it has a dramatic effect to 
the persistent currents of n x = 51, 52. In the second and third rows of Fig. HJ the 
sinusoidal-like variations on the left become sawtooth-like on the right, and have much 
larger magnitudes (similar to the system with n x = 50). This, of course, is due to the 
energy level crossing of the newly formed edge-modes. 

We now focus on the edge spin current in the presence of SOI. Spin-resolved edge 
currents In a (I La) are shown on the right (left) side of Fig. [5j The spin-resolved 
edge current has a discontinuity at $ = 0, similar to those in Fig. HI However, in 
Fig. [5l the sawtooth curves ride on background values of ±1 (i.e. opposite spins move 
along opposite directions). Therefore, it is a more stable transport with respect to the 
change of magnetic flux. As n y gets larger, even though the jump is smaller, the overall 
magnitudes of these currents remain close to one. 

A few symmetries can be observed in Fig. [51 For example, Il+ = Ir-, and 
II- = Ir+- Also, Ir+(4>) = — I r- (—(/)), similarly for the left edge. Notice that the 
charge current for the right edge, Ir = Ir + + Ir- , is small but nonzero due to the slight 
asymmetry between Ir±(—<P) and Ir±{4>). On the other hand, the spin current for the 
right edge, 1^ = Ir + — Ir-, is roughly of value 2 (in units of t/An). 

In the figure, the persistent spin current is nonzero even if there is no magnetic flux. 
Such a current, unlike persistent charge current, does not require the breaking of time- 
reversal symmetry. However, the edge states being helical is essential. Furthermore, 
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since the magnitude of the edge spin current is not reduced in a larger ring (in contrast 
to PCC), it is possible to have it in a large ring (if phase coherence around the ring could 
be preserved). For a flat sheet with open boundary (and no external bias), the edge spin 
current cannot persist. A continuous spin circulation without external driving force is 
possible only in a ring-like structure. In principle, the moving spins would induce an 
electric dipole [H] (might be oscillating in time) that can be felt. In practice, detecting 
such an induced electric field remains a great experimental challenge. 

3.2. Narrow ribbon 

As the ribbon becomes narrower, the edge states on the two sides will couple with each 
other [26]. Such a coupling may open the Dirac point (when n x ^ 3n + 2) at the Fermi 
energy. Fig. [6] shows the oscillatory energy gaps for different widths of the ribbon. One 
can see that the system remains semi-metallic when n x = 3n + 2, with or without SOI 
[27] . When n x ^ 3n + 2, as n x becomes smaller, inter-edge coupling produces a larger 
gap (triangles in Fig. |6]). For reference, the bulk gap at k y = is also plotted (solid 
squares). Notice that the bulk gap is larger at n x = 3n + 2, when the energy gap for 
the edge state vanishes. 

For a narrow ribbon, the edge modes are no longer well-defined. Therefore, we focus 
on the total persistent currents carried by the whole ring. Fig. [7J shows the persistent 
currents of spin-up and spin-down electrons in a narrow ring. As one can see, since 
I + — J_, the charge current I is twice of the values in Fig. [7J while the spin current 
I s is zero. There could be remnant spin currents along the left and the right edges, 
but they would cancel each other. For small widths (the left column), the currents for 
n x ^ 3n + 2 can barely be seen |28j. This is due to the finite energy gaps shown in 




Figure 6. Energy gaps for narrow ribbons with n y = 100. Solid circles are the energy 
gaps for A50 = 0. The system is semi-metallic for n x — 3n + 2, and is insulating 
otherwise. When A50 = 0.05, we show both the edge gaps (triangles) and the bulk 
gaps (solid squares). 



0.03 




Figure 7. Persistent currents for electrons in a narrow ring with n v = 100 and 
Aso — 0.05. The left column shows the persistent currents for up and down spins 
when n x = 14, 15, 16, while the right column is for n x = 29, 30, 31. 



Fig. |6j The case of n x = 3n + 2 is special since the Dirac point remains closed. Such a 
dramatic contrast is less apparent in the right column. 

3.3. Zeeman interaction 

The analysis so far applies to the cases when the electron in the ring feels the Aharonov- 
Bohm effect through the vector potential of the magnetic flux, but not the magnetic 
field itself. In the following, we consider a more realistic situation when the graphene 
ring itself is immersed in a magnetic field. As a result, the electrons can couple with 
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Figure 8. The energy levels near k y = for a ring with n x — 50, n y — 100, and 
^so — 0.01 are plotted with respect to <j> near <f> = 1. Each curve for g = is separated 
to two curves when g is nonzero. 
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Figure 9. Persistent current for a graphene ring with n x — 50, n y = 100, and 
Xso — 0.01. Insets show magnified views near the jumps. 



the field through their spins. Some order-of-magnitude estimate is in order: For a ring 
with n y -sites, the cross- sectional area is (3/47r)n^a 2 (a = 2.456A). In order to have one 
flux quantum inside the ring, the magnetic field needs to be B ~ 2.87 x 10 5 /n 2 (in units 
of T). The Zeeman gap at one flux quantum is 33.25/n 2 eV (g = 2). In comparison, 
the energy difference between neighboring fc^-states in the Brillouin zone (see Fig. [2]) 
is roughly of the order of 2t/n y (t ~ 2.8 eV). Therefore, the Zeeman gap is relatively 
unimportant, unless n y is of the order of ten or less. 

Fig. |8] shows the energy levels near the energy gap (with and without the Zeeman 
interaction). When g = 0, there are two parabolas with a minute energy gap due to the 
inter-edge coupling. When g ^ 0, the upper (and lower) parabola separates into two 



Persistent currents in a graphene ring with armchair edges 12 




^R 

Figure 10. Reduction of the edge spin currents (in units of t/An) due to the Rashba 
coupling. The parameters are n x = 50, n y — 100, Xso — 0.05, and the Zeeman 
coupling is not considered. 

parabolas. The separation is proportional to the strength of the Zeeman interaction. 
As a result, it is possible for the upper and lower parabolas to intersect with each other 
and produce two crossings (two Dirac points). This will alter the qualitative behavior 
of the 7(0) curve. 

The Zeeman effect on the PCC can be seen in Fig. [9j Near an integer flux quantum, 
a vertical jump splits to two (inset), due to the two Dirac points in Fig. [HJ Such a split 
is proportional to the magnitude of the magnetic field. As a result, the I ((f)) curve is no 
longer periodic in the magnetic flux. This is the main effect of the Zeeman interaction. 

3.4- Rashba interaction 

For a graphene sheet, the Rashba coupling is possible only if the mirror symmetry with 
respect to the plane is broken, for example, by an electric field or a substrate. For a 
graphene tube, the Rashba coupling could be induced, for example, by an electric field 
transverse to the axis of the tube. Kane and Mele found that, if X R > 2y / 3Aso, then 
the graphene would no longer be a QSH insulator [2]. For a smaller Xr, even though 
the graphene is still in the QSH phase with helical edge states, the electron spins in 
the edge channels are no longer conserved. Therefore, the magnitude of the edge spin 
current is reduced. It is interesting to investigate such a reduction in the spin transport. 
The ring geometry proves to be most convenient since one does not have to worry about 
the disturbance from external leads. 

Fig.lTOlshows the magnitude of the edge spin currents. When Xr = 0, the magnitude 
is very close to 2. It gradually decreases to zero at the critical value of Xr = 2y / 3Aso ) 
beyond which the graphene is semi- metallic [2] and there is no more edge spin current. 
Overall, the spin currents for spin-x and spin-?/ (calculated from semiclassical expressions 
similar to Eq. OH])) remain small for the whole range of the Rashba coupling. One can 
see that, even if Xr is as large as Xso, the spin current is still dominated by the spin- z 
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component. The slight deviation from the ideal values of 2 and on both ends is a 
result of the finite-size effect and inter-edge coupling. Similar reduction of the spin in 
edge states, when the spin is not conserved, has been observed in various spin-resolved 
measurements of Angle-Resolved Photo-Emission Spectroscopy on three-dimensional 
topological insulators [29] . 

4. Conclusion 

By adding a spin-orbit interaction to a graphene ribbon, one can generate helical edge 
states along armchair edges. When the ribbon is wrapped to form a ring that encloses a 
magnetic flux, there are persistent charge and spin currents. We studied the persistent 
currents in rings with different radii and widths. The bulk state contribution to the 
currents vanishes if the radius of the ring is large, but the edge state contribution 
persists. For broad ribbons, the decay length of the edge mode at the Dirac point 
is inversely proportional to Xso- The persistent spin current for one helical edge is 
roughly t/2n. For narrow ribbons, inter-edge coupling between the edge states opens 
the Dirac point (when n x ^ 3n + 2) and reduces the persistent currents. For most of 
the parameters being studied, the Zeeman effect is too small to have much influence. 
On the other hand, the Rashba coupling is found to deteriorate the edge spin transport, 
which vanishes beyond a critical value when the graphene is no longer a QSH insulator. 
Our work shows that the tubular QSH insulator driven by a magnetic flux is an ideal 
configuration for theoretically exploring the edge transport. 

Finally, a remark on the effect of electron-electron interaction, which is beyond 
the scope of this work. In Ref.[23], it has been shown that the Hubbard type, short- 
ranged interaction could induce edge magnetization. This is related to the following 
fact: the honeycomb lattice is bi-partite, and the Hubbard interaction tends to induce 
antiferromagnetic order. For a zigzag edge, since the lattice points along the edge are 
from the same sub-lattice, their spins would point to the same direction. That is the 
main reason why a zigzag edge is ferromagnetic. This mechanism will not work well 
for the armchair edge as it contains lattice points from both sub-lattices. Therefore, 
magnetization induced by the Hubbard interaction is less likely for the armchair edge. 
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